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Single-photon which is initially uncorrelated with atom, will evolve to be entangled with the atom 
on their continuous kinetic variables in the process of resonant scattering. We find the relations 
between the entanglement and their physical control parameters, which indicates that high entan- 
glement can be reached by broadening the scale of the atomic wave or squeezing the linewidth of 
the incident single-photon pulse. 

PACS numbers: 03.65.Ud, 42.50.Vk, 32.80.Lg 



I. INTRODUCTION 



Quantum entanglement is of fundamental importance 
in the theory of quantum nonlocalitypj as well as 
in quantum information 2j. Recently, photon-atom 
entanglement is frequently discussed in their finite 
Hilbert spacesQ, such as, the polarizations of photon 
or the internal states of atom. With the progress of 
micro-cavity quantum electrodynamics Q and high 
coupling artificial atompj, single photon raises its ability 
to affect considerably not only the atom's internal state 
but also its external motion. As a result, it gives rise 
to some basic questions related to the photon-atom 
entanglement on their infinite kinetic degree of freedom. 

In recent studies 0)0, entanglement in the continuous 
kinetic variables between single-photon and atom is 
mostly discussed in the process of single-photon emission 
with atomic recoil, where the atom is initially pumped 
to its excited level and the single-photon is prepared 
"intrinsically" by the atomic spontaneous emission. In 
our work, however, the resonant single-photon is initially 
injected from a tuneable single-photon generator [3, 
whereas an artificial atom is placed freely in vacuum on 
its steady state ("artificial" indicates that the atomic 
coupling to the single-photon is stronger than usual, 
which ensures the interaction observable 9]). We find 
that, after the interaction, the scattered single-photon 
will be entangled to the atom at a higher degree 
compared with the case of solely spontaneous emission. 
We explain this phenomena as the coherent pumping 
of the incident photon and evaluate it with a defined 
"entanglement pumping coefficient" . 



*Author to whom correspondence should be addressed. E-mail: 
hongguo@pku.edu.cn, phone: +86-10-6275-7035, Fax: +86-10- 
6275-3208. 



To describe the degree of entanglement, firstly, we use 
the ratio (R) between the conditional and unconditional 
variance in momentum to evaluate the two particles' 
correlation in the probability amplitude of their wave 
function, which is experimentally accessible and can be 
seen as the "amplitude entanglement" in momentum 
space ^JEHi secondly, we use the standard Schmidt 
decomposition and treat Schmidt number k\u$ as 
a criterion for the full entanglement contained both 
in amplitude and phase. For both criterions R and 
K, we revealed their dependencies on the physical 
control parameters r and 77, and compare them in some 
region of interests, from which it is shown that: higher 
entanglement can be achieved by either broadening the 
scale of the atomic wave or squeezing the linewidth of 
the incident single-photon. Transmitted photon is also 
considered, which is different to the scattered photon, 
and exhibits little entanglement with the atom due to its 
interference with the transparent wave (initially incident 
photon wave profile). 



II. THEORETICAL ANALYSIS 

As shown in Fig. 1 (a), the two-level atom with 
transition frequency uj a and mass m is placed freely 
in vacuum, the ground and excited states of which 
are denoted by |1) and |2), respectively. The incident 
single-photon from some generator is resonant with the 
atom and exhibits a superposed state of different fock 
states due to its linewidth. For the realistic consideration 
in some experiments^), we fi x the photon detector and 
atom detector in opposite directions and make them 
both in the x-z plane for simplicity as in Fig. I (b), 
the angle 6 can be chosen to observe the scattering in 
needed directions. 

Under the rotating wave approximation (RWA) the 
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FIG. 1: (a) Single-photon interacts resonantly with free two- 
level atom. 

(b) The incident photon is scattered by the atom, angle 9 is 
fixed to determine the direction of the detection. 

(c) Schematic diagram for the absorption-emission process. 
The process of emission with atomic recoil will generate en- 
tanglement between the recoiled atom and the scattered pho- 
ton due to the momentum conservation. 



Hamiltonian can be written in Schrodinger picture as: 
{hp) 2 

fc 

+ hJ2 [gik^ale-^ + H.c.] , (1) 
fc 

where hp and f denote atomic center-of-mass mo- 
mentum and position operators, denotes the atomic 

operator — 1, 2), ag and at are the annihilation 

and creation operators for the light mode with photonic 
wave vector fc and frequency — ck, respectively. Note 
the summation is performed over all coupled modes 
in the continuous Hilbert space. We also suppress 
the polarization index in the summation as well as in 
photon state, since we can always choose a particular 
polarization to detect the photon. g(k) is the dipole 
coupling coefficient. 

As there is only one photon in the interaction, the basis 
of the Hilbert space can be denoted as \q, lg, i) (i = 1, 2), 
where the arguments in the kets denote, respectively, the 
wave vector of the atom, and of the photon, and the 
atomic internal state. At time t the state vector can 
therefore be expanded as: 

|V) = ^d(g,fc, +^C 2 (g,t)|g,0,2). (2) 

q,k Q 

Substituting Eqs. (1) and (2) into Schrodinger equa- 



tion yields: 

iA(q, k, t) = g(k)B(q + fc, ty[ck-^a-^{2q+k)-k\t^ ( 3 ) 
iB(q, t)=J2 9*0)A(q- k, k, tyK-cfc+^(2?-fc)-fc> 

k 

(4) 

where A, B are the slowly varying parts of C\ and C2, 
i.e.: 

A(q,k,t) = C 1 (q,k,t)e i ^ +Ck ^ t , (5) 
B(q,t) = C^K^+^X. (6) 

Suppose the atom is initially in the ground state and 
has zero average velocity, the initial condition can be set 
as: 

A(q, k,t = 0) = XoG(q)P(k - k ), (7) 
B(q,t = 0)=0, (8) 

where G(q) = G x (q x )G y (q v )G z (q z ) and P(k - k ) = 
Px{k x )Py{ky)P z {k z — k ). In this case, functions Gi(qi) 
and Pi{ki) (i = x,y,z) have zero center value and 
bandwidths Sqi and Ski separately. The coordinates are 
chosen as in Fig. 1 (b), where we make the incident 
direction as z-axis. xo is the normalized factor and 
fco = (0,0, ^) is the resonant wave vector. 

We proceed to solve the equations with Laplace trans- 
formation and single pole approximation 15j and yield: 

B(q,t) = -i X0 J29*(k) x (9) 

G(q- k)P{k~k ) ^a-ck+J^(2q^k)-k]t _ e -iU-Tt\ 

iL + T + i[uj a -ck+^{2q-k)-k} 

where the frequency shift L and atomic linewidth T are 
given as: 

L = y \g(k)\ 2 

r = 7rJ2\9(k)\ 2 S(u a -ck). 

k 

We can simplify Eq. (9), by replacing the term G(q — k) 
with G(q— kg) since the momentum bandwidth Sqi due 
to the recoil is normally much larger than the photon 
linewidth Ski] also, we can replace k with ko in the term 
2^-(2<f — fc) • fc. With these approximations, the first 
term in the curly bracket can be seen as the antifourier 
transform of the product of the photonic shape and 
Lorentzian shape, and will cause a decay at a time scale 
max{jr,-i-}; the second decay term e - lLt ~ Tt i s due 
to the spontaneous emission. Then, one can directly 
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find that B(q,t — ► oo) — > 0. In the further calculations, 
we ignore the frequency shift since it can be treated 
as a modification of the atomic transition frequency, 
and regard the slowly varying function g(k) as a constant. 



With the approximations mentioned above, from Eqs. 
(3) and (9), we obtain the steady solution of A(q, k,t — > 
oo): 



A(q, k, t -> co) = X oG(q)P(k - fc ) + 



E 



Xo|<?| 2 G(q + fc-/c ) 



c£;-w a - [h(q + k)} 2 /2mh+(hq) 2 /2mh 
P(h - fc ) 



ck — ck\ 



JL-(2q + k)-k + ±(q + k)-k -^kt 



(10) 



From Eq. (10), one sees that the final state is a superpo- 
sition of the transparent wave (initially incident photon 
wave profile, depicted by the first term on the r.h.s.) 
and scattering wave (second term on the r.h.s.). In the 
scattering part, the atom and the photon are entangled 
due to the process of photon absorption and emission 
with atomic recoil, which is sketched in Fig. 1 (c). One 
may find that the Lorentzian-Gaussian factor in the 
scattering part is very similar to that in the case of 
spontaneous emission with recoil|6|, where the Gaussian 
term is a reflection of momentum conservation and the 
Lorentzian term indicates the energy conservation. 

The general formula (10) can be used to analyze the 
photon scattered in different directions. Without loss 
of physical generality, we choose the initial conditions 
for the atom as Gi(qi) — e~( qi > Sqi ^ , and for the pho- 
ton P,(fcj) = 1/ (ki/Ski + 1) which is exactly the case if 
the incident single-photon is generated by spontaneous 
emission. As a remark, we point out that all the conclu- 
sions in the following keep available when the incident 
photon is chosen to be other shapes such as Gaussian or 
whatever. 



(a) (b) 
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Ak z Ak z 

FIG. 2: (a) and (b) are contour and density plots of |j4^| 2 
with the condition r z = 1, rj x = 10; (c) and (d) are contour 
and density plots of \Aq\ 2 with the condition r z — 1, r/ z — 10. 



III. AMPLITUDE ENTANGLEMENT IN 
SCATTERED PHOTON 

To make the physical results more evident and avoid 
unnecessary mathematical complexity, we focus our at- 
tention on the photon scattered perpendicular to the in- 
cident direction, i.e., = \. Then we project Eq. (10) 
into the subspace \(q x ,0,0)) ® |l(fc x ,o,Q))! with the same 
approximations used in Eq. (9), and yield: 

, _ N-exp [-(Ag x -^Ak x ) 2 / v l] 

{Ak x + Aq x + ^ + i) [(Ak x + Aq x )/r z + i] 
„ iV-cxp^Afe/Tfa) 2 ] 

(Ak x + Aq x + i) [{Ak x + Aq x )/r z + i] ' 



where Ak % = A Ql = - k ), Vi = 

Ti = (i = x,y,z) are all defined dimensionless 

parameters. Note that rj x and t z contain all the 
physical parameters that determine the nature of the 
atom-photon system, thus can be treated as physical 
control parameters for the atom and the photon, re- 
spectively. We neglect tiny terms in Eq. (11) due to 
fifcg <C mT and hk <C mc in realistic conditions. N is 
the normalization factor where iV 2 = y/2(l + t z ) / T z r\ x . 

From Eq. (11) and Fig. 2, one sees that, variables Aq x 
and Ak x play the symmetric role in the two Lorentzian 
functions. It makes the probability amplitude |^4|-| 2 lo- 
calized along the diagonal of the momentum space, which 
implies the nonfactorization of the photon-atom wave 
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FIG. 3: (a) Relation between R and the two control parame- 
ters (r z , rj x ). (b) Sectional views of (a), with r) x = 1, 5, 10, 20 
from bottom to top. The ratio R is calculated from variable 
Aq x with Ak x fixed at the origin. 



function, and then will generate entanglement between 
the two particles. In fact, we can treat the ratio (R) 
of the conditional and unconditional variances for Aq x 
or Ak x as an evaluation of entanglement 01 ■ This ra- 
tio, compared to the Schmidt number K, reveals more 
obvious analytic dependence for the entanglement on its 
control parameters r\ x and t z , and is also experimentally 
directly accessible jl6| . 

We proceed to calculate the ratio for variable Aq x , i.e., 
R = 5Aq x nElc /SAq x omc , where the unconditional vari- 
ance is obtained from the single-particle observation as: 



(5 2 A(? si„ gl c = {Aq 2 } _ {Aqx) 2 

= J dAk x dAq x Aq 2 x \Az 



(12) 



dAk x dAq x Aq x \Ai 



and coincidence measurement gives the conditional vari- 
ance at some specified Ak x : 



5 2 Aq°°™ = (A<£)a*. 

_ $ dAq x Aql\A^ 



- (Aq x 
( J dAq x Aq x \A f \ 2 " 

"I /dA&|A 5 P 



(13) 



Substituting Eqs. (11)-(13) into the definition of R, we 
yield R{i] x , t z ) as a function of parameters rj x and t z , the 
result of which is illustrated in Fig. 3 with Ak x fixed 
at the origin. From that, one can see that the entan- 
glement increases monotonously when r\ x increases or r 2 
decreases, which indicates that higher entanglement can 
be achieved by squeezing the linewidth of the incident 
photon or broadening the wave packet of the atom. In 
particular, when r\ x > 1, we have: 



R 



(14) 



from which it is found that the entanglement increases 
linearly with r\ x and will be abruptly enhanced when r z 
tends to zero. As a remark, we emphasize that all the 



conclusions above hold qualitatively the same either if 
Ak x is specified otherwise or one calculate the ratio R 
from the other variable Afc^. 

The ratio i?, which can be obtained experimentally 
by comparing the momentum dispersion variance, is 
an appropriate quantification for the entanglement 
contained in the probability amplitude correlation (thus 
can be seen as an evaluation of the "amplitude entan- 
glement"). Next, we can see that it reveals a correct 
varying tendency for the entanglement with its control 
parameters. However, the definition of R is dependent 
on its representation space and different choices for 
the basis of Hilbert space will cause distinct values of 
R. This is because we only use the amplitude of the 
wavefunction to construct i?, and then all entanglements 
included in phase ^l| is lost. 

To obtain the "total entanglement", we calculate the 
Schmidt number fl3| and compare it with the entangle- 
ment ratio R in the following section. 



IV. FULL ENTANGLEMENT IN SCATTERED 
PHOTON 



Mathematically for a bipartite system in pure state, 
the entanglement of an unfactorable wavefunction can be 
completely characterized by the Schmidt number, which 
is denoted by K = E^Lo^n),— ; wnere ^'n s are eigen- 
values of the integral equation [12||: 

JdAk' xP p (Ak x , Ak x )<f>JAk x ) = \ n <f> n (Ak x ), (15) 

the density matrix for photon is defined as: 

P p {Ak Xl Ak' x ) ee fdAq x A f (Aq x , Ak x )A%(Aq x , Ak' x ), 

(16) 

where, note that we have taken away the time-dependent 
phase in the density matrix since it does not contribute 
to entanglement. Although we do it with the photon, 
Schmidt number can be equally obtained through the 
atomic density matrix, and the cigenfunctions of atom 
[ifj n (Aq x )] can be related to those of photon through: 



^n(Aq x ) 



1 



dAk x A^{Aq x ,Ak x )(t>* n {Ak x ), (17) 



where <j) n {Ak x ) and ^(Ag^.) (n — 1. 2 ■ •■) form complete 
orthonormal sets for the photon and atom respectively. 
With these discrete modes, the unfactorable wavefunc- 
tion can be expanded into a sum of factored products 
uniquely: 



.4 



E 



<Kip n (Aq x )(t> n {Ak x ). (18) 
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FIG. 4: Schmidt number if and the amplitude entanglement 
degree R in dependence on r z with r) x = 10. Spots are numer- 
ical results for K whereas solid line is plotted for R. The inset 
shows them as functions of r/ x with r z fixed, lines from bottom 
to top are depicted as: R(t z = 10), K(t z = 10), R(t z = 1), 
K{t z = 1), J?(r z = 0.1), K(t z = 0.1), respectively. 



Then, the Schmidt number K, which is an estimation 
of the number of modes that are "important" in making 
up the expansion of Eq. (18), serves as a quantitive 
measurement of entanglement [7] 13] . Note K is in- 
dependent from representation since all A's keep the 
same in different representations, thus can be seen as 
a quantity of the full entanglement information (both 
amplitude and phase entanglement) kept in the collective 
wavefunction. 

Since Eq. (15) is not analytically solvable, we use a 
discrete eigenvalue equation to approximate the integral 
equation. Up to a reliable precision, we use 1000 x 1000 
matrices to carry out the diagonalization, and collect 
some of the results in Fig. 4, where we also compare 
Schmidt number K with the amplitude entanglement 
ratio R. 
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FIG. 5: Entanglement pumping coefficient EPC as a function 
of t z . The solid line is its fitted function 1.1/t* + 1.5. 



(f] ^S> 1). This difference indicates that, although in both 
cases, entanglement is generated from momentum conser- 
vation in the process of photon emission with atomic re- 
coil, the absorption of the incident photon will add some 
entanglement due to its coherent pumping effect. As K is 
linear with rj (or rj x ), we define the "entanglement pump- 
ing coefficient" as: 



EPC 



slope of K{rj x ) in scattering 
slope of K (jj) in spontaneous emission 



since the constant term in K(rf) plays a minor role 
when entanglement is large. The defined coefficient 
EPC shows the times that entanglement is increased 
by the coherent pumping of an incident photon. As it 
is independent on the atomic parameter, it reflects the 
ability of entanglement of the photon separately. We 
collect some numerical results in Fig. 5 and fit it with 
EPC « l.l/r z + 1.5 within r z S (0, 1), from which, one 
sees that EPC increases rapidly when r z diminishes, 
which also implies that, if the incident photon is pre- 
pared monochromatically on its limit condition, i.e., 
r 2 — > 0, the scattered photon will be highly entangled to 
the recoiled atom. 



From the numerical results, we find that, similar to 
the ratio R, K rises linearly with parameter r\ x and 
will increase rapidly when the linewidth of incident 
photon is squeezed narrower to the atomic linewidth 
r, i.e, t z < 1; secondly, when r z is fixed, the slope of 
K(Vx) is always larger than that of R(rj x ), which means 
that more entanglement information will transfer to 
phase when rj x becomes larger, and this phenomena will 
become more evident when r z is reduced, e.g., when 
t z = 0.1, R « 1.587k + 1-39 whereas K m 3.44t7 x + 0.08, 
which indicates that more than half of the entanglement 
information will be unavailable by momentum dispersion 
observation when rj x goes large on this condition. 

Another phenomena is notable, when t z = 1, i.e., the 
linewidth of the incident photon is not squeezed and can 
be prepared directly by spontaneous emission from the 
same atom, we find K ss 0.7577a, + 0-16 (Vx 3> 1) whereas 
in the case of spontaneous emission[(| K w 0.28?7 + 0.72 



We plot the amplitude of the first three Schmidt 
modes for the photon with r\ x = 10 and r z = 1 in Fig. 
6. We find that their number of peaks in momentum 
space is proportional to the Schmidt mode index, but 
the separations of different peaks are more distinct than 
in the case of spontaneous emission|7j. 



V. TRANSMITTED PHOTON 

To consider the transmitted photon, we make the ob- 
servation angle 9 = 0, and yield the collective wavefunc- 
tion from Eq. (10): 

A Q = -XoG z (q z )P z (k z ) (19) 
n ( T V r x T y G z (q z )P z (k z ) 
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FIG. 6: First three Schmidt modes for the scattered and 
transmitted photon. Left column is for the scattered photon 
with t z = 1 and r) x = 10; right column is for the transmit- 
ted photon with t z = 1, r) x = 10, and j(-^) 2 T x r y = 1 for 
illustration. 

One can see that, in Eq. (19), the first term describes 
that the two particles are free of interaction and keep 
their initial factorable wave form; the second term 
reflects the entanglement. Usually, the second term is 
much smaller than the first one since (^) 2 "C 1, but 
one can enlarge it by choosing some special physical 
system, such as the artificial atom with low excited 
level and high coupling to its resonant modes. However, 
this improvement can add few entanglement between 
the transmitted photon and recoiled atom, because 
interference between the two terms in Eq. (19) will 
weaken the correlation of the two particles at a great 
deal. To make it clear, we show the contour and density 
plots for the probability amplitude of Aq in Fig. 2 
on an artificial condition \{-^) 2 T x T y — 1, and yield 
R w K < 2 in this situation. 

The eigenfunctions of transmitted photon for the first 
three modes with r\ z = 10 and t z — 1 are collected in Fig. 
6, from which one can see that, due to the interference, 
the corresponding modes of the transmitted photon ex- 



hibit one peak less than that of the scattered photon. 



VI. CONCLUSION 

We analyze the physically fundamental interaction 
between a single photon and a free artificial atom in 
vacuum. With a few physical approximations, the 
general solution of the photon-atom wave function 
is obtained, from which, it is found that the initially 
uncorrelated particles will evolve to be entangled due 
to momentum conservation in scattering. To evaluate 
the entanglement in the scattering, firstly, we use an 
experimentally accessible parameter R, which denotes 
the ratio between momentum variance in single-particle 
and in coincidence observations, and yield its simple 
dependences on the two physical control parameters 
f] x = and t z = secondly, we use standard 

Schmidt decomposition to reveal the full entanglement 
information and find out its varying tendency similar 
to that of R, which indicates that high entanglement 
can be achieved by cither squeezing the linewidth 
of the incident photon or broadening the scale of 
atomic wave packet. Furthermore, compared with 
spontaneous emission, we defined a parameter EPC 
to evaluate the entanglement enhancement due to the 
coherent pumping effect of the resonant incident photon. 
In the end, we found out that, for the transmitted 
photon, one can expect little entanglement due to the 
interference between the transparent and scattered wave. 
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I. INTRODUCTION 



Quantum entanglement is of fundamental importance 
in the theory of quantum nonlocalitypj as well as 
in quantum information 2j. Recently, photon-atom 
entanglement is frequently discussed in their finite 
Hilbert spacesQ, such as, the polarizations of photon 
or the internal states of atom. With the progress of 
micro-cavity quantum electrodynamics Q and high 
coupling artificial atompj, single photon raises its ability 
to affect considerably not only the atom's internal state 
but also its external motion. As a result, it gives rise 
to some basic questions related to the photon-atom 
entanglement on their infinite kinetic degree of freedom. 

In recent studies 0)0, entanglement in the continuous 
kinetic variables between single-photon and atom is 
mostly discussed in the process of single-photon emission 
with atomic recoil, where the atom is initially pumped 
to its excited level and the single-photon is prepared 
"intrinsically" by the atomic spontaneous emission. In 
our work, however, the resonant single-photon is initially 
injected from a tuneable single-photon generator 
whereas an artificial atom is placed freely in vacuum on 
its steady state ("artificial" indicates that the atomic 
coupling to the single-photon is stronger than usual, 
which ensures the interaction observable). We find that, 
after the interaction, the scattered single-photon will be 
entangled to the atom at a higher degree compared with 
the case of solely spontaneous emission. We explain 
this phenomena as the coherent pumping of the incident 
photon and evaluate it with a defined "entanglement 
pumping coefficient" . 



To describe the degree of entanglement, firstly, we use 
the ratio (R) between the conditional and unconditional 
variance in momentum to evaluate the two particles' 
correlation in the probability amplitude of their wave 
function, which is experimentally accessible and can be 
seen as the "amplitude entanglement" in momentum 
space ^JEHi secondly, we use the standard Schmidt 
decomposition and treat Schmidt number k\u$ as 
a criterion for the full entanglement contained both 
in amplitude and phase. For both criterions R and 
K, we revealed their dependencies on the physical 
control parameters r and ry, and compare them in some 
region of interests, from which it is shown that: higher 
entanglement can be achieved by either broadening the 
scale of the atomic wave or squeezing the linewidth of 
the incident single-photon. Transmitted photon is also 
considered, which is different to the scattered photon, 
and exhibits little entanglement with the atom due to 
its interference with the transparent wave. 



*Author to whom correspondence should be addressed. E-mail: 
hongguo@pku.edu.cn, phone: +86-10-6275-7035, Fax: +86-10- 
6275-3208. 



II. THEORETICAL ANALYSIS 

As shown in Fig. 1 (a), the two-level atom with 
transition frequency uj a and mass m is placed freely 
in vacuum, the ground and excited states of which 
are denoted by |1) and |2), respectively. The incident 
single-photon from some generator is resonant with the 
atom and exhibits a superposed state of different fock 
states due to its linewidth. For the realistic consideration 
in some experiments^), we S x the photon detector and 
atom detector in opposite directions and make them 
both in the x—z plane for simplicity as in Fig. 1 (b), 
the angle 9 can be chosen to observe the scattering in 
needed directions. 

Under the rotating wave approximation (RWA) the 
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FIG. 1: (a) Single-photon interacts resonantly with free two- 
level atom. 

(b) The incident photon is scattered by the atom, angle 9 is 
fixed to determine the direction of the detection. 

(c) Schematic diagram for the absorption-emission process. 
The process of emission with atomic recoil will generate en- 
tanglement between the recoiled atom and the scattered pho- 
ton due to the momentum conservation. 



Hamiltonian can be written in Schrodinger picture as: 
{hp) 2 

fc 

+ hJ2 [gik^ale-^ + H.c.] , (1) 
fc* 

where hp and f denote atomic center-of-mass mo- 
mentum and position operators, denotes the atomic 

operator — 1, 2), ag and at are the annihilation 

and creation operators for the light mode with photonic 
wave vector fc and frequency — ck, respectively. Note 
the summation is performed over all coupled modes 
in the continuous Hilbert space. We also suppress 
the polarization index in the summation as well as in 
photon state, since we can always choose a particular 
polarization to detect the photon. g(k) is the dipole 
coupling coefficient. 

As there is only one photon in the interaction, the basis 
of the Hilbert space can be denoted as \q, lg, i) (i = 1, 2), 
where the arguments in the kets denote, respectively, the 
wave vector of the atom, and of the photon, and the 
atomic internal state. At time t the state vector can 
therefore be expanded as: 

|V) = ^d(g,fc, +^C 2 (g,t)|g,0,2). (2) 

q,k Q 

Substituting Eqs. (1) and (2) into Schrodinger equa- 



tion yields: 

iA(q, k, t) = g(k)B(q + fc, t yWk-^a-Jt(2q+k)-k]t^ (3) 
iB(q, t) = J2 9*0)A(q- k, k, tyK-cfc+^(2?-fc)-fc> 

k 

(4) 

where A, B are the slowly varying parts of C\ and C2, 
i.e.: 

A{q,k,t) = C 1 (g ) fc,t)e i ^+ Cfe )*, (5) 
B(q,t) = CMty^+^t. (6) 

Suppose the atom is initially in the ground state and 
has zero average velocity, the initial condition can be set 
as: 

A(q, k,t = 0) = xoG(q)P(k - k ), (7) 
B(q,t = 0)=0, (8) 

where G(q) = G x (q x )G y (q y )G z {q z ) and P(k - k ) = 
Px{k x )Py{ky)P z {k z — k ). In this case, functions Gi{qi) 
and Pi{ki) (i — x.y,z) have zero center value and 
bandwidths Sqi and Ski separately. The coordinates are 
chosen as in Fig. 1 (b), where we make the incident 
direction as z-axis. xo is the normalized factor and 
fco = (0,0, — ) is the resonant wave vector. 

We proceed to solve the equations with Laplace trans- 
formation and single pole approximation 15j and yield: 

B(q,t) = -i X oY,9*(k) x (9) 

fc 

G(q- k)P(k-k ) je'K-^+^^-fe)-^ _ e -iLt-rt\ 
iL + T + i[u a - ck+ ^{2q-k) ■ k] 

where the frequency shift L and atomic linewidth T are 
given as: 

^ uj a ~ck + ^(2q-k)-k ' 
T = n^2\g(k)\ 2 S(u a -ck). 

k 

We can simplify Eq. (9), by replacing the term G(q — k) 
with G(q — fco) since the momentum bandwidth Sqi due 
to the recoil is normally much larger than the photon 
linewidth Ski ; also, we can replace fc with fco in the term 
2^7 (2<f — fc) • fc, for the very sharp decay resonance of 
the denominator. With these approximations, the first 
term in the curly bracket can be seen as an antifourier 
transformation on the production of the photonic shape 
and Lorentzian shape, and will cause a decay at a time 
scale max{^,^-}; the second decay term e - li *- r * j s 
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due to the spontaneous emission. Then, one can directly 
find that B(q,t — + oo) — > 0. In the further calculations, 
we ignore the frequency shift since it can be treated 
as a modification of the atomic transition frequency, 
and regard the slowly varying function g(k) as a constant. 



With the approximations mentioned above, from Eqs. 
(3) and (9), we obtain the steady solution of A(q, k,t — > 
oo): 



A(q, k, t -> oo) = X oG(q)P(k - k Q ) + 



E 



xo\g\ 2 G(q + k - k ) 



ck-oJa- [h(q + k)} 2 /2mh+{hq) 2 /2mh 
P{ki - fco) 



^ z [cfc - cfel _ £.( 2 q+ fc) • fc + ±(q + fc) ■ fc - ^fc 2 



(10) 



From Eq. (10), one sees that the final state is a superpo- 
sition of the transparent wave (first term on the r.h.s.) 
and scattering wave (second term on the r.h.s.). In the 
scattering part, the atom and the photon are entangled 
due to the process of photon absorption and emission 
with atomic recoil, which is sketched in Fig. 1 (c). One 
may find that the Lorentzian-Gaussian factor in the 
scattering part is very similar to that in the case of 
spontaneous emission with recoil|6|, where the Gaussian 
term is a reflection of momentum conservation and the 
Lorentzian term indicates the energy conservation. 

The general formula (10) can be used to analyze the 
photon scattered in different directions. Without loss 
of physical generality, we choose the initial conditions 
for the atom as Gi(qi) — e~( qi > Sqi ^ , and for the pho- 
ton Pi(ki) = 1/ (ki/Ski + 1) which is exactly the case if 
the incident single-photon is generated by spontaneous 
emission. As a remark, we point out that all the conclu- 
sions in the following keep available when the incident 
photon is chosen to be other shapes such as Gaussian or 
whatever. 



(a) (b) 
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lk z Ak z 

FIG. 2: (a) and (b) are contour and density plots of | 2 
with the condition r z = 1, rj x = 10; (c) and (d) are contour 
and density plots of \Aq\ 2 with the condition r z = 1, r\ z — 10. 



III. AMPLITUDE ENTANGLEMENT IN 
SCATTERED PHOTON 

To make the physical results more evident and avoid 
unnecessary mathematical complexity, we focus our at- 
tention on the photon scattered perpendicular to the in- 
cident direction, i.e., = \. Then we project Eq. (10) 
into the subspace \(q x ,0,0)) ® |l(fc x ,o,Q))! with the same 
approximations used in Eq. (9), and yield: 

, _ N-exp [-(Ag x -^Ak x ) 2 / v l] 

{Ak x + Aq x + ^ + i) [(Ak x + Aq x )/r z + i] 
„ iV-cxp^Afe/Tfa) 2 ] 

(Ak x + Aq x + i) [{Ak x + Aq x )/r z + i] ' 



where Ak % = A Ql = ^fe - fc ), Vi = 

Ti = (i = x,y,z) are all defined dimensionless 

parameters. Note that rj x and t z contain all the 
physical parameters that determine the nature of the 
atom-photon system, thus can be treated as physical 
control parameters for the atom and the photon, re- 
spectively. We neglect tiny terms in Eq. (11) due to 
fifcg -C mT and hkg -C mc in realistic conditions. N is 
the normalization factor where iV 2 = ^/2(l + T z )/7r2T z r] x . 

From Eq. (11) and Fig. 2, one sees that, variables Aq x 
and Ak x play the symmetric role in the two Lorentzian 
functions. It makes the probability amplitude |^4|-| 2 lo- 
calized along the diagonal of the momentum space, which 
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FIG. 3: (a) Relation between R and the two control parame- 
ters (r z , rj x ). (b) Sectional views of (a), with r) x = 1, 5, 10, 20 
from bottom to top. The ratio R is calculated from variable 
Aq x with Ak x fixed at the origin. 



implies the nonfactorization of the photon-atom wave 
function, and then will generate entanglement between 
the two particles. In fact, we can treat the ratio (R) 
of the conditional and unconditional variances for Aq x 
or Ak x as an evaluation of entanglement 01 ■ This ra- 
tio, compared to the Schmidt number K, reveals more 
obvious analytic dependence for the entanglement on its 
control parameter s ru and r z , and is also experimentally 
directly accessible fl6| . 

We proceed to calculate the ratio for variable Aq Xl i.e., 
R = 5 Aq x nElc J '5 Aq x olnc , where the unconditional vari- 
ance is obtained from the single-particle observation as: 



(5 2 Ag singlo 



(Ad) - <A<fe) 2 

dAk x dAq x Aq 2 x \A« | 2 

dAk x dAq x Aq x \A i \ 



(12) 



and coincidence measurement gives the conditional vari- 
ance at some specified Ak x : 



(5 2 A(? coinc = (A(? 2 )a ^ _ {Aq: 



z/Afc. 



(13) 



JdAq x Aql\A; 

JdAq x \A^ 



' J dAq x Aq x \Ai 
JdAq x \AJ 



Substituting Eqs. (11)— (13) into the definition of R, we 
yield R{i] x , t z ) as a function of parameters rj x and t z , the 
result of which is illustrated in Fig. 3 with Ak x fixed 
at the origin. From that, one can see that the entan- 
glement increases monotonously when r\ x increases or r 2 
decreases, which indicates that higher entanglement can 
be achieved by squeezing the linewidth of the incident 
photon or broadening the wave packet of the atom. In 
particular, when r\ x > 1, we have: 



R 



Vx 



2^ 



(14) 



from which it is found that the entanglement increases 
linearly with r\ x and will be abruptly enhanced when r 2 



tends to zero. As a remark, we emphasize that all the 
conclusions above hold qualitatively the same either if 
Ak x is specified otherwise or one calculate the ratio R 
from the other variable Ak x . 

The ratio R, which can be obtained experimentally 
by comparing the momentum dispersion variance, is 
an appropriate quantification for the entanglement 
contained in the probability amplitude correlation (thus 
can be seen as an evaluation of the "amplitude entan- 
glement"). Next, we can see that it reveals a correct 
varying tendency for the entanglement with its control 
parameters. However, the definition of R is dependent 
on its representation space and different choices for 
the basis of Hilbert space will cause distinct values of 
R. This is because we only use the amplitude of the 
wavefunction to construct R, and then all entanglements 
included in phase is lost. 

To obtain the "total entanglement", we calculate the 
Schmidt number [Tflj and compare it with the entangle- 
ment ratio R in the following section. 



IV. FULL ENTANGLEMENT IN SCATTERED 
PHOTON 

Mathematically, for a bipartite system in pure state, 
the entanglement of an unfactorable wavefunction can be 
completely characterized by the Schmidt number, which 
is denoted by K = (Yl™=otfi)ZLi wnere A„s are eigen- 
values of the integral equation jl2| : 



J dAk' xP p (Ak x ,Ak' x )M^K 



X n <p n (Ak x ), (15) 



the density matrix for photon is defined as: 

P p (Ak x ,Ak' x ) = J dAq x A i (Aq x ,Ak x )Al(Aq x ,Ak' x ), 

(16) 

where, note that we have taken away the time-dependent 
phase in the density matrix since it does not contribute 
to entanglement. Although we do it with the photon, 
Schmidt number can be equally obtained through the 
atomic density matrix, and the cigenfunctions of atom 
[ip n (Aq x )] can be related to those of photon through: 



ip n {Aq x ) 



1 



dAk x A«(Aq x ,Ak x )ct>* n (Ak x ), (17) 



where (j) n (Ak x ) and tp n (Aq x ) (n = 1, 2 • ••) form complete 
orthonormal sets for the photon and atom respectively. 
With these discrete modes, the unfactorable wavefunc- 
tion can be expanded into a sum of factored products 
uniquely: 

-4§(A 9x ,Afc x ) = ^2 \ZKip n (Aq x )(j) n (Ak x ). (18) 
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K(t z ) & R(r z > 
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FIG. 4: Schmidt number if and the amplitude entanglement 
degree R in dependence on r z with r) x = 10. Spots are numer- 
ical results for K whereas solid line is plotted for R. The inset 
shows them as functions of r/ x with r z fixed, lines from bottom 
to top are depicted as: R(t z = 10), K(t z = 10), R(t z = 1), 
K{t z = 1), J?(r z = 0.1), K(t z = 0.1), respectively. 



Then, the Schmidt number K, which is an estimation 
of the number of modes that are "important" in making 
up the expansion of Eq. (18), serves as a quantitive 
measurement of entanglement [7] 13] . Note K is in- 
dependent from representation since all A's keep the 
same in different representations, thus can be seen as 
a quantity of the full entanglement information (both 
amplitude and phase entanglement) kept in the collective 
wavefunction. 

Since Eq. (15) is not analytically solvable, we use a 
discrete eigenvalue equation to approximate the integral 
equation. Up to a reliable precision, we use 1000 x 1000 
matrices to carry out the diagonalization, and collect 
some of the results in Fig. 4, where we also compare 
Schmidt number K with the amplitude entanglement 
ratio R. 



a. 
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FIG. 5: Entanglement pumping coefficient EPC as a function 
of t z . The solid line is its fitted function 1.1/t* + 1.5. 



(f] ^S> 1). This difference indicates that, although in both 
cases, entanglement is generated from momentum conser- 
vation in the process of photon emission with atomic re- 
coil, the absorption of the incident photon will add some 
entanglement due to its coherent pumping effect. As K is 
linear with rj (or rj x ), we define the "entanglement pump- 
ing coefficient" as: 



EPC 



slope of K{rj x ) in scattering 
slope of K (jj) in spontaneous emission 



since the constant term in K(rf) plays a minor role 
when entanglement is large. The defined coefficient 
EPC shows the times that entanglement is increased 
by the coherent pumping of an incident photon. As it 
is independent on the atomic parameter, it reflects the 
ability of entanglement of the photon separately. We 
collect some numerical results in Fig. 5 and fit it with 
EPC « l.l/r z + 1.5 within r z S (0, 1), from which, one 
sees that EPC increases rapidly when r z diminishes, 
which also implies that, if the incident photon is pre- 
pared monochromatically on its limit condition, i.e., 
r 2 — > 0, the scattered photon will be highly entangled to 
the recoiled atom. 



From the numerical results, we find that, similar to 
the ratio R, K rises linearly with parameter r\ x and 
will increase rapidly when the linewidth of incident 
photon is squeezed narrower to the atomic linewidth 
r, i.e, t z < 1; secondly, when r z is fixed, the slope of 
K(Vx) is always larger than that of R(rj x ), which means 
that more entanglement information will transfer to 
phase when rj x becomes larger, and this phenomena will 
become more evident when r z is reduced, e.g., when 
t z = 0.1, R « 1.587k + 1-39 whereas K m 3.44t7 x + 0.08, 
which indicates that more than half of the entanglement 
information will be unavailable by momentum dispersion 
observation when rj x goes large on this condition. 

Another phenomena is notable, when t z = 1, i.e., the 
linewidth of the incident photon is not squeezed and can 
be prepared directly by spontaneous emission from the 
same atom, we find K ss 0.7577a, + 0-16 (Vx 3> 1) whereas 
in the case of spontaneous emission[(| K w 0.28?7 + 0.72 



We plot the amplitude of the first three Schmidt 
modes for the photon with r\ x = 10 and r z = 1 in Fig. 
6. We find that their number of peaks in momentum 
space is proportional to the Schmidt mode index, but 
the separations of different peaks are more distinct than 
in the case of spontaneous emission|7j. 



V. TRANSMITTED PHOTON 

To consider the transmitted photon, we make the ob- 
servation angle 9 = 0, and yield the collective wavefunc- 
tion from Eq. (10): 

A Q = -XoG z (q z )P z (k z ) (19) 
n ( T V r x T y G z (q z )P z (k z ) 
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FIG. 6: First three Schmidt modes for the scattered and 
transmitted photon. Left column is for the scattered photon 
with t z = 1 and r) x = 10; right column is for the transmit- 
ted photon with t z = 1, r) x = 10, and j(-^) 2 T x r y = 1 for 
illustration. 

One can see that, in Eq. (19), the first term describes 
that the two particles are free of interaction and keep 
their initial factorable wave form; the second term 
reflects the entanglement. Usually, the second term is 
much smaller than the first one since (^) 2 "C 1, but 
one can enlarge it by choosing some special physical 
system, such as the artificial atom with low excited 
level and high coupling to its resonant modes. However, 
this improvement can add few entanglement between 
the transmitted photon and recoiled atom, because 
interference between the two terms in Eq. (19) will 
weaken the correlation of the two particles at a great 
deal. To make it clear, we show the contour and density 
plots for the probability amplitude of Aq in Fig. 2 
on an artificial condition \{-^) 2 T x T y — 1, and yield 
R w K < 2 in this situation. 

The eigenfunctions of transmitted photon for the first 
three modes with r\ z = 10 and t z — 1 are collected in Fig. 
6, from which one can see that, due to the interference, 
the corresponding modes of the transmitted photon ex- 



hibit one peak less than that of the scattered photon. 



VI. CONCLUSION 

We analyze the physically fundamental interaction 
between a single photon and a free artificial atom in 
vacuum. With a few physical approximations, the 
general solution of the photon-atom wave function 
is obtained, from which, it is found that the initially 
uncorrelated particles will evolve to be entangled due 
to momentum conservation in scattering. To evaluate 
the entanglement in the scattering, firstly, we use an 
experimentally accessible parameter R, which denotes 
the ratio between momentum variance in single-particle 
and in coincidence observations, and yield its simple 
dependences on the two physical control parameters 
f] x = and t z = secondly, we use standard 

Schmidt decomposition to reveal the full entanglement 
information and find out its varying tendency similar 
to that of R, which indicates that high entanglement 
can be achieved by cither squeezing the linewidth 
of the incident photon or broadening the scale of 
atomic wave packet. Furthermore, compared with 
spontaneous emission, we defined a parameter EPC 
to evaluate the entanglement enhancement due to the 
coherent pumping effect of the resonant incident photon. 
In the end, we found out that, for the transmitted 
photon, one can expect little entanglement due to the 
interference between the transparent and scattered wave. 
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